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ABSTRACT 



We have extended our earlier work on space weathering of the youngest S-complex 
asteroid families to include results from asteroid clusters with ages < 10 6 years and to 
newly identified asteroid pairs with ages < 5 x 10 5 years. We have identified three S- 
complex asteroid clusters amongst the set of clusters with ages in the range 10 5-6 years 
- (1270) Datura, (21509) Lucascavin and (16598) 1992 YC2. The average color of the 
objects in these clust ers agrees with the color predicted by the space weathering model of 



Willman et al. I ( 120081 ). SDSS 5-filter photometry of the members of the very young asteroid 
pairs with ages < 10 5 years was used to determine their taxonomic classification. Their types 
are consistent with the background population near each object. The average color o f the S - 
complex pairs is PC\ = 0.49 ±0.03, over 5a redder than predicted by lWillman et al. I ( 120081 ). 
This may indicate that the most likely pair formation mechanism is a gentle separation due 
to YORP spin-up leaving much of the aged and reddened surface undisturbed. If this is the 
case then our color measurement allows us to set an upper limit of ~ 64% on the amount of 
surface disturbed in the separation process. Using pre-existing color data and our new results 
for the youngest S-complex asteroid clusters we have extended our space weather model to 
explicitly include the effects of regolith gardening and fit separate weathering and gardening 
characteristic timescales of t w = 960 ± 160 My and r g = 2000 ± 290 My respectively. The 
first principal component color for fresh S-complex material is PC\ = 0.37 ± 0.01 while the 
maximum amount of local reddening is APCi = 0.33 ±0.06. Our first-ever determination of 
the gardening time is in stark contrast to our calculated gardening time of r g ~ 270 My based 
on main belt impact rates and reasonable assumptions about crater and ejecta blanket sizes. 
A possible resolution for the discrepancy is through a 'honeycomb' mechanism in which the 
surface regolith structure absorbs small impactors without producing significant ejecta. This 
mechanism could also account for the paucity of small craters on (433) Eros. 



Keywords: Asteroids, surfaces; Spectrophotometry; Spectroscopy 



-4- 



Introduction 



After three decades of development our understanding of the spa ce weathering phe- 
nome non may be converging on its cause and effects. Originally proposed (j Chapman and Salisbury 
19731 ) as a solution to the mismatch between the color of the most common meteorites, or- 
dinary chondrites, and their likely source region, inner main belt S-complex asteroids, the 
space weathering hypot hesis is that the su rface colors of asteroids change with exposure to 
the space environment (jHapkd 1 19731 . Il968l ) . The idea has held up to the test of time and in 
recent years has led to measurements of the rate of color change on S-co mplex asteroids yield 



ing empirical models that predict the rate o f reddening of their surfaces ( jJedicke et al. 



2004 



Nesvorny et al. Il2005t IWillman et al. 11200a ). This work extends our space weathering model 
to the youngest dated asteroids in the main belt that are less than 10 6 years old and, for 
the first time, measures the surface gardening rate on the asteroids due to impact-generated 
regolith cycling. 



Willman et al. I (120081 ) 's space weathering model applies exclusively to S-complex aster- 



oids. Other types of asteroids may also undergo space weathering and we expect that they 
would obey a model of similar functional form with different parameters. For instance, if C- 
complex asteroids undergo space weathering they must have a different color range because 
they occupy a different region of principal co mponent space and th e weathering rate and 
even the sense of coloration may be different (INesvorny et al. II2005I ). The gardening rate 
would probably also be different due to different material density and strength. 

In this work we study the space weathering of S-complex asteroids without particular 
concern for the agent causing the weathering though the likely culprit is energy deposition due 



to par ticle bombardment e.g. micrometeorites, solar protons and cosmic rays. iMarchi et al 



(120061 ) claim that Sun-related effects are dominant but they assumed that micrometeorite 
impacts are independent o f heliocentric d istance. Since neither their flux nor their velocity 
are distance independent (jCinatalal 119921 ) it is unclear whether it is justified to downgrade 
their contribution to the space weathering agent inventory. Laboratory-based studies suggest 
that the asteroidal space weathering mechanism, long confused with the lunar processes pro- 
ducing agglutinates, involves the coating of near surface s emi-transparent nanophase grains 
by a vapor or sputter deposited film of metallic iron {e.g. 

2oooh . 



Sasaki et al. 2001; Pieters et al 



While measuring the color of asteroid surfaces is relatively straightforward and asteroid 
taxonomy based on colors has been a mainstay of planetary science for decades the determi- 
nation of asteroid surface ages is a relatively recent innovation. The surface ages of asteroid 
family members may be determined from dynamical simulations of the evolution of the fam 



ily's orbit element distribution (e.g. IVokrouhlicky and Nesvornyll2006a ; IVokrouhlicky et al. . 
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2006bl ; INesvorny et al. 1120051 . 120021 ; iMarzari et al. 1 1 1 9 9 5h . The dynamical methods include 
family size frequency distribution (SFD) modeling, global main belt SFD modeling, mod- 
eling of family spreading via thermal forces and backward numerical integration of orbits. 
The combination of asteroid surface colors with their ages allowed the first astron omical 
determination of the space weathering rate (IJedicke et al. Il2004t iParker et al. 1 120081 1 . 



Space weathering rate measurements utilize the fact that asteroid family members are 
fragments of a single collisionally disrupted parent body that formed at the same time and 
assume that fresh unweathered material from the collision debris cloud coats all members of 
a family with a homogeneous regolith layer. Therefore, all family members should start out 
with similar color and weather in tandem. 

We ignore grain size color effects for a number of reasons. First, we believe that unless 
grain size is correlated with age its effect will average out over the members of a family. Sec- 
ond, we will see below that km scale asteroids, in contrast to moon sized bodies with highly 
comminuted regoliths, are gravitationally sorted rubble piles whose surfaces are dominated 
by boulders. Fines are largely absent, perhaps sequestered beneath the bouldered surface, 
so asteroid regolith appears to vary with size. This would diminish grain size observ a tiona l 
effe cts in the sub-millige e environment of small asteroids. Finally, INesvorny et al. I (120051 ) 
and I Jedicke et al. I (120041 ) did not find any correlation between color and size, just color and 
age. 



Jedicke et al. I (120041 ) 's space weathering rate measurements encapsulated 'the color' of 



an asteroid's surface as its first principal component color 
PC l = 0.396 (u-g) + 0.553 (g - r) + 0.567 (g 



+ 0.465 (g - z) 



T 



that correlates strongly with the average slope of the spectrum. The principal component 
color provides the linear combination of filter magnitudes having the greatest variability over 
the sample of asteroids. In this case PC\ is specific to the asteroids in the Sloan Digital Sky 
Survey (SDSS) 2 nd Moving Object Catalog. While most of the color variability in the SDSS 
asteroid sample is captured in PCi, the second principal component, 



PC 2 = -0.819 (u -g) + 0.017 (g 



0.090 (g 



0.567 (g 



(2) 



( INesvorny et al. 1 120051 ) corresponds to the asteroid spectrum's curvature and is used in this 
analysis to identify asteroid taxonomy. 

Motivated by the fact that uniform s urface irradiati o n sho uld produce an exponentially 
decaying amount of unweathered surface IJedicke et al. I (J2004J) showed that the color of S- 
complex asteroids as a function of time t can be expressed as 



Pd(t) = PCi(0) + APCi[l - e 



-{t/rY 



(3) 
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where PC i(0) is the unweathered color of fresh surface material, APC± is the magnitude 
of the weathering color change after a long period of time, r is the charac teristic time for 
space weathering and a is a generalizing exponent. IWillman et al. I ( 120081 ) improved upon 
the earlier work by refining the color of the < 5 My old Iannini family, eliminating the Eos 
family that is no longer considered to be in the S-complex, and carefully refitting the data 
to the functional form given above to determine PC±(0) = 0.31 ±0.04, APC\ = 0.31 ±0.07, 
r = 570 ± 220 My and a = 0.53 ± 0.19. 

The characteristic time scale of 570 ± 220 My for color change in th e main asteroid belt 
( IWillman et al. 1 120081 ) is in agreement with pulsed laser experiments ( iSasaki et al. I l200lh 
on silicate pellets intended to simulate micrometeorite bombardment. They suggested a 
characteristic time for weathering of 100 My at 1 AU, equivalent to about 700 My in the main 
belt assuming that the Sun is the source of the weather ing agent and a r~ 2 dependence on 
its strength. However, ( iLoeffler et al. 1120091 ) argue that ISasaki et al. I (120011 ) overestimated 
the space weathering characteristic timescale due to an incorrect flux calculation leading. 
Basically, the time scale extrapolation from lab results is complicated because it is unclear 
how to evaluate the contribution from several different factors. For instance, in the main 
belt micrometeorite impact speeds will be lower (~ 5 km/s) than at 1 AU (~ 20 km/sec) 
and their impact energy will be ~16x weaker. 



Pieters et al. I (120001 ) estimated characteristic aging times of 100 — 800 My for lunar 



surfaces by comparing ages from craters dated radiometrically or by cosmic ray exposure 
ages to spectral differences. Correcting the rate to the center of the main belt suggests 
weathering times on the order of ~600-4800 My. Even though the lunar surface is not 
identical to S-complex asteroid surfaces the time scales are similar. 

Fin ally, we note that cra ters on (243) Ida are bluer than their surrounding background 



terrain (Veverka et al. 



regolith (IVeverka et al 



99$) . The craters correspond to freshly exposed and unweathered 



while other parts of the asteroid's surface indicates an age of 



about 1 Gy ( iGreenberg et al. Ill996l ). The wide range in diameters (a proxy for crater age) 
of blueish crater suggests that the space weathering time must be long. 

On the o ther hand there are claims of faster surface weathering time scales such as 



Takatol ( 120081 ) 's upper limit of 450 kyr based on a shallow 1 fim absorption band observed 
on (1270) Datura. However, we expect that the space weathering phenomenon is a rel- 
atively subtle effect easily masked by stochastic variations between asteroids due to i.e. 
mineralogical and morphological differences and/or collisional and cratering events. The 
space weathering effect can only be identified in specific regions on an asteroid as observed 
with in situ spacecra ft measurements on (95 1) Gaspra (|Helfenstein et al. 1119941 ). (243) Ida 



( IVeverka et al. Ill996[ ). and (25143) Itokawa (llshiguro et al 



20071 )] or as an ensemble effect 
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on a statistically large sample of asteroids. It is therefore difficult to make a general con- 
clusion based on a single asteroid such as (1270) Datura. Indeed, our meas urement of that 
astero id's PC\ = 0.41 ± 0.02 is redder by 2.2<r than its predicted color from lWillman et al 



( 120081 ). Takato's measurement of a shallow 1 /im absorption band which would tend to red- 
den the overall spectrum is therefore generally consistent with our measurement. We expect 
that the 'redness' of individual small rubble pile asteroids in a sample is affected more by 
random surface variations than larger regolith-rich asteroids with finer surface materials. In 
this case (1270) Datura is redder than the Datura cluster average so its particular value can 
be misleading. The four (1270) Datura members in this work have mean PC\ = 0.305 with 
a RMS of 0.278. The large standard deviation is not due to measurement error — it is a 
result of intrinsic color differences between the members of the (1270) Datura family and 
is typical of other families. Therefore we would discount the resulting 450 kyr weathering 
time. 

Another sho rt time scale measure ment was reported from recent lab experiments on 
olivine powder by lLoeffler et al. I (120091 ) (and references therein) simulating solar wind effects 
at 1 AU by bombardment with 4 keV protons. They conclude that spectral reddening caused 
by the solar wind should saturate in ~ 5 kyr. 



(jVernazza et al. 1120091 ) combined spectral data from a sample of four members of young 
clusters with archival meteorite spectra and found that space weathering is substantially 
complete in < 1 My but then co ntinues at a slow e r pac e up to several Gy. They attribute 
the first stage to the solar wind (jStrazzulla et al. I ( 120081 ) lab simulation) and the second to 
micrometeorite bombardment. This scenario depends critically on the starting color that 
they derive from meteorite data. 

It is not clear how to reconcile the information from these various lab results that differ 
in time scale by several orders of magnitude with our observations. 



Willman et al. I ( 120081 ) 's weathering model was derived from families that were sev- 
eral My (e.g. Karin, Iannini) to several Gy (e.g. Eunomia, Maria) old. The family age 
estimates h ave typical errors of ~ 40% resulting from fundamental limitations in the dating 
techniques (jNesvorny et al. 1120071 ). 



Although there were nine known S-complex families with ages from tens of Mys to a 
few Gys at that time there were only two younger than 10 My. Refining the space weath- 
ering rate at even younger ages requires a large sample of young asteroid families, which 
are typically produced by the catastrophic disruption of small asteroids, have only a small 
number of detectable members, and are therefore difficult to identify. But in recent years 
the number of catalogued asteroid orbits has reached hundreds of thousands and includes 
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many asteroids smaller than 1 km. This large sample of asteroids al lows the identific a tion o f 
rare small cluster Jj] originat i ng; fro m collisions less than 1 My ago. Nesvorny et al. ( 2006 ); 
Nesvorny and Vokrouhlickyl (120061 ) found four such clusters with a total of 16 members. The 
clusters are named after their largest known members: (1270) Datura, (14627) Emilkowalski, 
(21509) Lucascavin and (16598) 1992 YC2. 



The progression to ever smaller clusters reached its logical limit with lVokrouhlicky and Nesvorny 

(120081 ) 's discovery of 60 pairs of asteroids in extremely similar orbits — much more alike than 



would be expected based on the density of proper elements for other asteroids with similar 
orbit elements. The pairs are thought to have formed less than 500 kyr ago based on the dy- 
namical evolut ion time of the pair member's orbits. Their list was updated and restricted to 
only 36 pairs (IPravec et al. II2009I ). Pairs belonging to known young families such as (1270) 
Datura or (832) Karin were excluded because of the possibility that the catastrophic impact 
and subsequent family formation process may create paired asteroids. Such cases were elim- 
inated to focus on pairs that were created in isolation and presumably by the same method. 
The formation method of pairs may involve a critical distinction from that of clusters or 
families. Whereas the latter form in catastrophic collisions, pairs have alternative possible 
formation methods. Some of these methods may be gentle processes not involving resetting 
of the surface as discussed in §6.31 



Shortly after the publication of IWillman et al. I ( 120081 ) we obtained observations of some 
members of the sub-My old families and found that they were significantly redder than pre- 
dicted. Did this imply a problem with our space weathering model or is the space weathering 
phenomenon more complicated than suggested by the simple model of eq. 13]? 

We were a l so con cerned with the functional form of the space weathering model of 
Willman et al. I (120081) i.e. eq. El I f spac e weathering is an isolated process the exponent a 
should be unity but IJedicke et al. I (120041 ) fit a = 0.53 ± 0.19. Their argument was that a is 
a generalizing factor that accounts for both space weathering and regolith gardening which 
acts to counteract the surface aging by slowing turning over the asteroid's surface. 

To investigate these questions we collected color and spectral data from several sources 
for members of t he young families and pairs. We obtained spectra of mem bers of the sub-My 
asteroid clusters ( INesvorny et al. 112006k INesvorny and Vokrouhlickyll2006l ) and were provided 



1 As pairs of asteroids arc (tautologically) composed of two asteroids, the established families have dozens 
to thousands of members, and known 'clusters' have between three and seven members we define a cluster 
as a small family with three to about ten asteroids. Pairs have distinct formation mechanisms from families 
and clusters. Families form from larger parent bodies than clusters and arc therefore older on average. The 
youngest known family, Iannini, is ^3 My old while the clusters are less than 1 My old. 
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spectra for some of the objects observed by iMothe-Diniz and Nesvornyl (120081 ). We also 
located archived photometry for 19 pair members in the Sloan Digital Sky S urvey Data 



Release 7 Moving Object Catalog 4 (SDSS DR7 MOC4) (IParker et al. 1 120081 ). We then 



developed a new asteroid surface color-age model that explicitly separates the effects of 
weathering and gardening and eliminates the need for the unexplained generalizing factor (a). 
The new model allowed us to measure the characteristic timescales for both weathering and 
gardening on the S-complex asteroids. Finally, we independently calculated the gardening 
timescale on main belt asteroids from their size distribution, impact rates and estimates of 
crater and ejecta blanket size. 



2. Space weathering versus regolith gardening effects 



The space weathering model of iJedicke et al. I (120041 ) summarized in eq. [3] uses a char 



acteristic time, r, along with a generalizing and unphysical exponent, a, to capture the 
time dependence of the gradual reddening of S-complex asteroid surfaces. The exponent was 
introduced because it was understood that there are more effects in play on an asteroid's 
surface than a single weathering component, e.g. multiple sources of space weathering with 
different time scales such as solar protons and ultraviolet radiation, micrometeorite bombard- 
ment and cosmic rays, in addition to the effect of regolith gardening which will counteract 
the space weathering. 'Gardening' is in part due to meteoroids that regularly strike the 
asteroids's surface and lift fresh sub-surface material to the regolith's top layer. Gardening 
may also take place through seismic shaking with subsequent regolith distribution or it may 
result from larger asteroid strikes which spread ejecta blankets beyond a crater. We do not 
distinguish between possible causes but lump them all under the banner of 'gardening' in the 
same way that all possible causes of color chang e in minerals are inc luded in 'weathering'. 



Here we extend the space weathering model of IJedicke et al. I (l2004f ) to explicitly include 
both weathering and gardening. 

Consider the relationship of unweathered surface, U, and weathered surface, W — by 
construction U + W = 1. Space weathering causes U to be replaced by W at the rate 

du- = _m (4) 

dt r w 

where r w is the characteristic space weathering time. Regolith gardening causes U to increase 
at the same rate that W decreases, 

dU + dW _ W(t) 

~df = ~dT ~ Tg ^ ' 
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where r g is the characteristic gardening time. 

The total rate of change of unweathered surface is then 



dU dU + dU~ 1 (\ 1 

+ -— = --[- + — )U. (6) 



dt dt dt Tg \Tg T„ 



Since U(0) 



g V T 9 Tw ) — |— 

U(r g ,T w ,t) = Tg (7) 



which has the desired properties that U(t = 0) = 1 and as t — > oo, U — > (1 + t 2 -)" 1 — the 
amount of unweathered asteroid surface after a long period of time is related to the relative 
rates of space weathering and regolith gardening. 

Replacing the unweathered surface term in eq. [3] with our new generalized unweathered 
surface term in eq. [7] yields 

PC x {t) = PCM + APd[l - U(r w ,r g ,t)} (8) 

which includes the effects of both weathering and gardening and will allow the two charac- 
teristic times to be separately determined by fitting to the asteroid families' color and age 
data. 

Previous attempts to explore the relationship between weathe red and unweathere d sur- 



face have foundered on the fundamental equations. For instance, iGault et al. I ( 119741 ) over- 



looked the restorative effects of regolith gardening and assumed that the fraction of surface 
that is undisturbed is simply a decaying exponential in time. 

We return to the model parameters that best fit the color-age data in §6.31 That garden- 
ing time determined from age-color data is compared to the resurfacing time determined from 
impacts in §6.41 In the next section we calculate the resurfacing rate using the independent 
method of impactors. 



3. Calculating the resurfacing rate from impacts 

The asteroid resurfacing rate due to impacts, S, is the product of the frequency of 
impacts, v(D), and the ejecta covered surface area, A e (D) per impact, where both terms are 



Gault et al. 



(|l974T ) eq. 1 is analogous to our cq. 4. 
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a function of impactor diameter D. Initially we place no limit on diameter which could extend 
from micrometeoroids up to parent body shattering size. The constraint on size will enter 
the discussion below in the lower limit of the impact integral. We realize that the asteroid's 



surfa ce may also be indirectly affected by impact-induced seismic shaking (IRichardson et al 



20051 ) w hich would inc r ease t he resurfacing rate. Seismic shaking could even happen without 



impact. iBinzel et al. I (120101 ) find that near Earth objects (NEOs) with likely close Earth 
encounters within the past 500 kyr have bluer surfaces than their counterparts that avoid such 
encounters suggesting that seismic shaking induced in this manner is effective at resurfacing. 
However, the mechanism of resurfacing during close encounters with massive planets is not 
important for MB asteroids. 

On the other hand, since asteroids in the 1 — 10 km diameter range have escape velocities 
of ~1 — 10 m/ s it is possible that a sizable fraction of ejecta on those asteroids is lost to space 
which would have the effect of decreasing the resurfacing rate. 

The impact frequency is the product of the differential size frequency distribution, N(D), 
the impact cross section, 

D T + D 



c(D T , D) = ( — — ^ ) (9) 



and the intrinsic collision probability, Pj = 2.86 x 10 -18 km -2 yr _1 (IBottke et al. I Il994l ) 
- the main belt wide average collision probability that a single member of the impacting 
population will hit a unit area of the target body per unit time. Dt is the target's diameter 
and the usual factor of n in eq. [9] is implicitly included in Pj. The cross section is determined 
by the limiting distance between centers at which the edge of the impacting asteroid just 
glances the edge of the target asteroid: Dt + d . 

Combining the above and summing over all impactor sizes up to the maximum diameter 
(D max ) that will not collisionally disrupt the target asteroid yields the surface gardening rate 
on the target asteroid: 

= Pi J ™ N{D) c(D T , D) A e (D) dD . (10) 

The resurfacing rate per unit target area is then S(Dt)/ttD t . The lower limit on the integral 
is ~ because small dust-size particles get blown away by solar radiation pressure and 
' dust' up to about 1 c m diameter spirals into the Sun due to the P oynting- Robert son eff ect 



( iDermott et al. II2002I ). The size-frequency distribution used here (IBottke et al. Il2005bl ) is 
only specified down to ~ 1 m diameter but we will show later that the size range from 1 cm 
to 1 m does not affect our conclusions. 



Since gardening of already gardened surface has no effect we are only interested in the 
fraction of gardening that takes place on weathered surface. Using the notation of §2] the 
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rate of gardening on weathered surface is: 

^l = w^il (11) 

leading to an exponentially decreasing amount of weathered surface with characteristic time 

r 9 = ^ (12) 
9 S(D T ) V ; 

Before combining equations M — H21 into final form there are several additional factors 
to consider that are covered in the following sub-sections: §3.11 the diameter of the largest 
impactor, D max , that will not shatter the target asteroid, §3.21 the diameter of the smallest 
impactor, D m i n , that will produce a significant ejecta blanket, §3.31 the diameter, D c , of the 
crater produced by an impactor of diameter, D, and §3.41 the diameter, D e , of the ejecta 
blanket. We finally calculate the regolith gardening rate in §3.51 



3.1. Largest non-shattering impactor, D max 

The size of an impactor that will shatter a target is determined by the shattering specific 
impact energy, Q* s (Dt). This is different from, and smaller than, the catastrophic disruption 
specific impact energy, Q* D (D T ), that applies to shattering the target and dispersing the 
fragments with high enough energy that reassembly into another rubble pile is impossible. 
Since we are concerned with regolith gardening we require that some portion of the asteroid's 
surface remain intact and use Q* s (Dt) to determine the diameter of the largest non-shattering 
impactor. Higher energy could, at a minimum, reset the entire asteroid surface to one with 
no weathering history or fission the target body. In either case the new surface would start 
with no history and would be unidentifiable with the previous surface. 

Smaller, typically monolithic objects are relatively resistant to catastrophic disrup- 
tion and, counterintuitively, large gravitational aggregates (rubble piles) have high collision 
strength because of their ability to absorb energy through non-elastic compression. While 
there are a wide variety of both theoretical and empirical Q* s (Dt) functions the specific 
energy is typically at a minimum for objects with diameters in the range 0.1 — 10 km (e.g. 



Holsapple et al. 1 120021 ). 



It can be shown that the ratio of impactor to target diameters required to shatter small 
asteroids is 

Dmax _ (2Q* S PT\ l/3 

D T V v 2 pi 



(13) 



13 



where v is the collision speed, p represents density and the sub-scripts / and T represent 
the impactor's and target's quantities respectively. The impactor and target densities are 
assumed to be equal as the most likely scenario in the inner main belt where our sample is 
located involves collisions between two S-complex asteroids of any size. 



Using iHolsapple et al. I (120021 ) 's Q* s (Dt) and v = 5 km/s typical of main belt asteroid 



collision speeds we find that D r 
in our data. 



,:Dr ~1:40 for asteroids of a few kms diameter like those 



3.2. Smallest ejecta blanket-creating impactor, D r 



We will show below that the calculated gardening rate is strongly dependent on the size 
of the smallest impactor capable of creating an ejecta blanket or, more generally, affecting 
an area substantially larger than its own cross section. Basically, because there are a large 
number of small impactors their contribution to the gardening rate could be substantial if 
not limited in same way by the target asteroid's structure or composition. For instance, if an 
asteroid's surface is mostly covered by meter-scale boulders it might require ^>-meter scale 
impactors to generate a significant gardening signature. 

The microgravity environment on small asteroids ( < 1 km diameter) displa ys surprising 
phenomena as illustrat ed by images of (25143 ) Itokawa ([Miyamoto et al. 1120071 ) obtained by 
the Hayabusa mission (IFujiwara et al. 1 120061 ) . The asteroid's surface is imbricated in many 
areas — cobble-sized stones coat the surface in a fairly regular pattern and are oriented 
in a common direction in the manner of bricks. Fine materials are largely absent from 
the surface except where they are collected in l ow gravity potential p onds. The fines may 
have sifted downwards through the larger rocks (IMivamoto et al. 1l2007l) in addition to being 



preferentially dispersed into space during impacts (1 Chapman! 119781 ). This rubble pile of 
gravel appears to have been size sorted into an inside out hierarchy; the largest rocks being 
at the surface and fines hidden inside. 

We think that an imbricated surface may act akin to medieval chainmail effectively 
warding off blows and protecting the surface from impac t cratering. Evidence supporting 
this mechanism may be found in IMivamoto et al. I ( 120071 ) Figures 2 e,f where light colored 
ejecta from the adjacent Komaba crater appears to coat one side of scattered bricks but 
there is no continuous fine coating of the surface beyond the crater rim. This suggests that 
the coating was not broadcast over the landscape as a powder but rather that these rocks 
were thrown from the crater to their current position carrying the coating with them. This 
would be the result if the coated rocks had originally overlain the position of an impact 
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which caused them to be ejected with a dust coating on their undersides. The key point 
is that the dust itself was inhibited from being broadcast to resurface a wide area. We 
therefore postulate a 'chainmail' mechanism that results in resistance to penetration and 
ejecta creation by impacto rs smaller t han so me multiple of the brick size. The mechanism 
was termed 'armoring' by IChapmanl (120021 ) but we believe 'chainmail' is a slightly more 
descriptive term for the proposed surface quality. Armor brings to mind rigid continuous 
sheets rather than the flexible linked elements in chainmail. The 'chainmail' mechanism may 
have a significant influence on the effect of regolith gardening on small asteroids since small 
meteoroid impacts will be relatively inefficient at gardening compared to the ejecta blankets 
produced by larger impactors. 

With this in mind we all ow that impa ctors create craters according to accepted crater 
scaling laws for D > D min (jMeloshl Il989l ). For smaller impactors, those in a size range 
corresponding roughly to the size of the imbricated surface regolith, we will set the diameter 
of the affected area to D. But what is the appropriate value for D m i n l 

The only two asteroids that have been imaged at a surface resolution better than ~50 
m are (433) Eros and (25143) Itokawa. Both show im bricated surfaces with only small areas 
covered by fines. Figure 4b of lMiyamoto et al. I ( 120071 ) provides the cumulative size frequency 
distribution of gravel of diameter d on Itokawa with N oc d~ 2 ' 8 but there is a bump near d = 
2 m which is confirmed by visual inspection of the boulder strewn areas and a broad hump 
near d = 0.3 m. The smooth areas dominated b y finer mat e rial a re the rare exceptions on 
(25143) Itokawa. For (433) Eros, Figure 3 from IChapmanl (120051 ) shows a surface that is 
'bumpy' near the resolution limit of 1 m which we think implies that the surface elements 
are ~1 m in size. 



With an ignorant assumption that the impactor needs to be > 10 x more massive than 
the surface elements to be large enough to apply the crater scaling laws it would imply 
that D m i n ~ 2 m. For the purpose of this work we scale that figure up by another order 
of magnitude and take D m i n = 4 m as our nominal case. We will show in the discussion 
that using smaller D min only serves to increase the discrepancy between our measured and 
calculated gardening rates so that the choice of the nominal value is unimportant. 



3.3. Crater diameter, D c 



While impact craters on Earth typically have diameters D c ~ 10 x D (jde Pater et al 



200 lh this relationship should probably not be be used for the low gravity high-porosity 



cratering impacts taking place between asteroids in the main belt. In this sub-section we 
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show that we can not use standard crater scaling laws because they appear to over-estimate 
the D c /D ratio on asteroids and we will argue that it makes more sense to use a D c j D ratio 
that is only weakly D-dependent. 



The general empirical scaling lawcleq. 5.26b of lde Pater et al. I (120011 ) for crater diameter 
is (in mks units) 

D c = 1.8 pf 11 Pt 1/3 9t°- 22 E° k 22 sin 1/3 9 D 013 (14) 



where gx is the gravitational acceleration on the target's surface, E k is the kinectic impact 
energy, and 9 is the angle of impact from the local horizontal. Since the impact energy, E k , 
is expressed in terms of D, pj and the impact speed v then 



D c = 1.34 ( ^V /3 g^ 22 v 0A4 sin 1 / 3 9 D 0J9 . 
\p T J 



(15) 



Ignoring the gravitational focussing on the Eartlfl, the ratio of crater diameters on an asteroid 
and the Earth for the same size impactor is 



D c (asteroid) 
D c (Earth) 



Pa' ^9a 



0.22 /v x 0.44 

v E 



(16) 



But can this result be used to determine the ratio over the three orders of magnitude differ- 
ence between the size of the Earth and the asteroids in our sample? 



Using p = 2.7 g/cm 3 with q = 0.46 cm/s 2 for (433) Eros (IKorvcanskv fc Asphaual2004f) 
and p = 2.6 ± 0.5 g/cm 3 (IBelton et al. Ill995l ) with g = 1 cm/s 2 (IKorycansky &: Asphaug 
20041 ) for (243) Ida yields a simple arithmetic mean target density for the two asteroids of 
2.65 g/cm 3 and gravity of 0.73 c m/s 2 . Assum i ng an impactor densit y equal to the average S- 
complex asteroid of 2 .63 g cm -3 (jHiltonll2002l ; iFujiwara et al. II2006I ). a density of 2.7 g cm -3 



fjCollins et al. 1 120051 ) for the upper layer of the Ear th's continental cru st, typical impact 
speeds on the Earth and asteroid of v ~ 20 km/s fjCollins et al. 1 120051 ) and v — 5 km/s 



fjBottke et al. Ill994l ) respe ctively, we find that Z ) e (asteroid)/.D c (Earth) ~ 2.7. Then, since 



D c ~ 10 x D on the Earth (jde Pater et al. Il200ll ) the scaling law leads to an analogous ratio 
for two S-complex asteroids of D c ~ 27 x D. 



However, iBottke et al. I (j2006bl ) have used crater counting methods and an assumed 
impactor size distribution to estimate that D c ~ 12.5 x D on (433) Eros and (243) Ida. So 
a decrease in scale from the Earth's diameter to the 10s of km size of (243) Ida and (433) Eros 



3 We repeated the calculation with a similar equation from lSchmidt et al. I (1987) having slightly different 
exponents which apply to porous media. The result did not alter our conclusions. 

4 We will show in the discussion that the systematic errors are larger than the amplitude of this effect. 
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only increases the D C :D ratio from 10 to 12.5 in contrast to the scaling calculation in the last 
paragraph that suggested a ratio of ~ 27 x. Thus, it appears that the scaling relation of eq. 
IT51 fails to capture essential physics involved in crater formation on small asteroids although 
it is not surprising that an extrapolation from the Earth's size do wn to kilomete r scale 
fails. We cannot use the Ide Pater et al. I (120011 ) equation or even the ISchmidt et al. I (119871 ) 
equation specific to porous media to calculate crater size on our few km scale asteroids. 

We empirically selected a power law to fit the two data points, Earth with mean radius 
6371 km and D c /D = 10 and Eros & Ida with mean radius 16.3 km and D c /D = 12.5 
yielding 

n f J) x -0.03350 

t - 1373 id) ■ (17) 

The near zero exponent indicates the weak dependence of the crater diameter ratio on target 
size and allows us to justify extending the relationship to the few km diameter range of the 
asteroids typical in our sample. 



3.4. Ejecta blanket diameter, D e 

The final key to calculating the asteroid regolith gardening rate from impact rates and 
cratering effects is the area affected by the crater ejecta. Like all the other terms in the 
calculations, this effect is difficult to characterize because cratering on small asteroids is not 
as well studied as on planets and their satellites. Due to their weak gravity a p ortion of 
ejecta from typical impacts, particularly the fine material ( Nakamura et al. 1 1994), may be 
thrown entirely clear of the target resulting in reduced ejecta coverage (jChapmanl 120021 ) . 



The fraction of material eject e d fro m large craters is lower for higher porosity objects 
( jHousen et al. 1120031 ). iBritt et al. I (120061 ) estimate macroporosities for coherent objects such 
as the moon, S-class asteroids, and C-class asteroids to be < 2 — 3%, 15 — 20%, and 25 — 50% 
respectively. Hence we would expect porosities of the asteroids composing our sample of 
S-complex families to be some 7x higher than that of the Moon leading to smaller ejecta 
blankets on asteroids compared to the Moon. 



Meloshl ( 119891 ) describes continuous ejecta blankets as typically coating the lunar surface 



out to over two crater rim radii from the center. His eq. 6.3.1 gives the radius of continuous 
ejecta, R e , as a nearly linear function of the crater rim radius, R c , with R e = (2.3 ±0.5) R l c m& 
for 1.3 km< R c <436 km. This function roughly applies to Mercury, Mars, the Jovian 
satellites Ganymede and Callisto, and the Saturnian satellites Dione and Rhea. 



Bearing in mind the uncertain porosity and gravity effects we use R e = (2.3 ±0.5) R 



1.006 
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as our baseline estimate and now attempt to account for the distal rays. We imagine a 
torus immediately outside the continuous ejecta disc with an area equivalent to int egrating 



the pa tchy distal ray coverage. Using the lunar crater Timocharis (see Figure 6.2 in iMelosh 



( 119891 ) ) as a prototypical case we estimate that replacing the rays with such a torus increases 
the continuous ejecta disc of ~ 2.3 R c to about 2.7 R c yielding an area ~ 38% larger than 
the continuous ejecta blanket alone. Converting radius to diameter and combining D e = 
(2.7 ± 0.5) Dl m& and eq. E] yields D e = (37.7 ± 7.0) D ime D~ om3r or D e ~ 38 D. 

Thus, the covered by ejecta for D > D min is 

A e (D T ,D) = - A V\ = HOO D-° mr D 2m . (18) 



For D < D min we simply assume that the affected area is equal to the impactor's cross 

1 



section: A e (D T ,D) = | D 2 . 



3.5. Resurfacing time distribution 



Combining all the discussion in this section above we find that the characteristic gar- 
dening time from eq. [12] is given by 



T g (N{D),D T ,Q* s ,v,pT,pi) 



9.8 X 10 14 Z&.067 



fD^ ax (Q* s ,v, PT , Pl ) Ar(D) (&±2\ D 201 dD 



(19) 



Figure [T] provides our nominal estimate of the characteristic gard ening time using 



Q* S {D T ) from Iliolsapple et al. I J2002|). v = 5 km/s JBottke et al. I Il994l) or = p T =2.63 



g cm" 3 f lHiltonll2002t iFuiiwara et al. 1120061 ) and N(D) from lBottke et al. I (l2005bh . We find 
that the gardening time frame is a mild function of diameter, only 23% slower on small (1 
km) asteroids than on large (100 km). The calculated resurfacing time for Dt ~ 4 km, the 
typical diameter of asteroids in our study, is ~ 270 My. Note that asteroids with diameters 
<700 m have disruption lifetimes shorter than either the calculated gardening time or the 
weathering time implying that few of them would be fully weathered or g ardened. Since 
(25143 ) Itokawa is ~535 m long it is unlikely to be completely reddened and llshiguro et al. 
( 120071 ) provide evidence of some space weathering. Unfortunately, the three bvw filter bands 
used in situ on (25143) Itokawa are not comparable to the five SDSS ugriz filters and we 
cannot calculate PC\ to determine where (25143) Itokawa's color falls on our scale. 



As new surveys detect asteroids of smaller diameters we may be able to discover this 
transition size, further constraining the weathering/gardening model. However, the asteroids 
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Fig. 1.- 

Calculated gardening time (solid curve) from asteroid impact and cratering rates as a function of 
target diameter. The shaded region is our formal estimate of the statistical error on the gardening 
given the uncertainty on all the input parameters. The actual systematic error is much larger as 
discussed in the text. Our measured color-derived gardening time is the si ngle data point in th e 



upper left. The asteroid disruption lifetime is given by the dashed curve (jBottke et al. Il2005bl ) 



The histogram (right ordinate) shows the size distribution of all the asteroids used in this study 
assuming an albedo for the S-complex asteroids of 0.20. 



in our sample have lifetimes > 2 Gy so their record of space weathering and gardening should 
not be distorted by an early demise. 

The above result for the resurfacing time determined by impacts is compared to the 
gardening time determined via age-color data in §6.41 
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4. Data acquisition and reduction 



We observed or obtained spectra or photometry for ten asteroids within the four known 
sub-My clusters as shown in Table [TJ The data came from three s ources: 1) new spec - 



tra and photome try were acquired by our team usin g Keckll/ESI (ISheinis et al. 1 12002| ) 



UH2.2 m/SNIFS flLantz et al\\2W4\ and IRTF/Spe X flRavner et al. Il2003h . 2) spectra were 
also provided by iMothe-Diniz and Nesvornyl (120081 ). hereafter refer red to as MDN, an d 3) 



we identified photometry for some sources from SDSS DR7 MOC4 (IParker et al. II2008I ). 



Our spectroscopic reduction techniques for the new data followed genera lly accepted 
practi ces using standard IRAF and IDL procedures that are explained in detail in lWillman et al. 
(120081 ). Procedures applied to some of the spectrographic data that may not be standard 
include: 



• using 2-D arc lamp spectra to create a transform map to straighten the 2-D asteroid 
and analog spectra. 

• median-combining straightened dome flats into a column normalized master flat. i.e. 
the average of all pixels in any column (at the same wavelength) in the master dome 
flat was fixed at unity to correct for pixel-to-pixel variations in the quantum efficiency 
at all points on the CCD at the same wavelength. 

• creating a water band model by dividing solar analog spectra from the same star that 
show the greatest difference in water band amplitudes. This allowed us to enhance the 
water bands for that solar analog. Combining such cases from different stars produced 
the master water band spectrum that allowed better cancell ation of atmospheric water 
absorption bands redward of 800 nm (IBus and Binzelll2002l ). 



• dividing the master water band spectrum into the asteroid spectrum and allowing the 
strength of the master to vary in order to minimize the distortion due to water bands. 



JBus and Binzelll2002h 



binning the spectra into 10 nm bins in the manner of IBus and Binzell (120021 ) such that 
the realized resolution was in the range 44 < R < 92 over the wavelengths spanning 
440 - 920 nm. 



The spectral reduction procedures generally apply to UH2.2 m/SNIFS and IRTF/SpeX 
spectra although some steps were performed automatically by the SNIFS pip eline. The 



IRTF /SpeX data reduction is facilitated by an IDL based package called Spextool ( jCushing et al. 
20041 ) which includes preparation of calibration frames, processing and extraction of spectra 
from science frames, wavelength calibration and flux calibration. 
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Table [TJ Sub-My Asteroid Clusters 



Asteroid 


Observation (HST) 


Source 


Mag (H/V) 


Type 


Axis 


(1270) Datura 


2007 Oct 28 


UH/SNIFS 


12.5/17.1 


SI 


2.234 


(1270) Datura 


2008 Mar 2 


Keckll/ESI 


12.5/16.5 


SI 


2.234 


(1270) Datura 


2007 Nov 16 


MDN 


12.5/16.7 


Sk 


2.234 


(90265) 2003 CL 5 


2007 Mar 15 


MDN 


15.4/19.3 


Sq 


2.235 


(60151) 1999 UZ 6 


2007 Mar 15 


MDN 


16.1/20.5 


Sk 


2.235 


(60151) 1999 UZ 6 


2001 Mar 18 


SDSS 


16.1/20.1 


Sq/Sk 


2.235 


(203370) 2001 WY 35 


2007 Sep 2 


MDN 


17.6/21.5 


O/Q 


2.235 


2003 UD112 


2003 Sep 25 


SDSS 


18.4/20.0 


Sq 


2.233 


(14627) Emilkowalski 


2006 Oct 1 


UH/SNIFS 


13.1/16.6 


T 


2.598 


(14627) Emilkowalski 


2008 Apr 13 


UH/SNIFS 


13.1/17.3 


T 


2.598 


(14627) Emilkowalski 


2008 Apr 18 


IRTF/SpeX 


13.1/17.4 




2.598 


(14627) Emilkowalski 


2004 Jan 27 


SDSS 


13.1/18.3 


D 


2.598 


(16598) 1992 YC2 


2007 Aug 17 


MDN 


14.7/20.3 


Sq 


2.621 


(16598) 1992 YC2 


2000 Jan 1 


SDSS 


14.7/17.0 


s 


2.621 


(21509) Lucascavin 


2006 Aug 1 


UH/SNIFS 


15.0/19.2 


Sk 


2.281 


(21509) Lucascavin 


2008 Apr 18 


IRTF/SpeX 


15.0/18.1 




2.281 


(180255) 2003 VM 9 


2008 Mar 2 


Keckll/ESI 


17.0/19.7 


Sk 


2.280 


(209570) 2004 XL 40 


2007 Aug 20 


MDN 


17.1/20.4 


Sq 


2.281 



Table 1: 

Observations and basi c pro perties of members of four sub-My old asteroid clusters adapted from 



Nesvornv et al. I ([20061 ) and 



Nesvornv and Vokrouhlickvl (|2006l ). The four clusters are separated by 



table lines and named after the first object listed in each cluster. Dates of observation and the 
data source are listed in the second and third columns respectively. Absolute magn i tude (H 



is 



listed along with the apparent magnitude (V) on the dat e of observation from I JPD. (12009) . We 



have determined the taxonomic type based on the SMASS ([Bus and Binzell 12002 : iBus et al 



2002) 



classification system as described in the text. The tight clustering in semi-major axis within each 
cluster is a consequence of their family membership. 



We used the solar colors of iBlanton et al. I ( 120071 ) in o ur PC^ and P C? measurement 
and absolute magnitude corrections for u and z bands from lSDSsl ( 2006b ). 



To combine the spectral results with the photometric results described below in a con- 
sistent manner we required a method to calculate PC\ (eq. [[]) for the sp ectra. To do so 



we rel ied on the relationship between spectral slope s and PC\ derived by IWillman et al. 
( [20081 ) . They used a sample of 133 asteroids common to both SMASS and the SDSS 3 rd 
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Data Release to determine that 

Pd = (0.87 ± 0.02) x (s/fim' 1 ) + (0.082 ± 0.012) (20) 

where the SMASS slope s is determined over the wavelength range from 0.44 — 0.92 /im as 
the best fit line pivoting through the normalization point at 0.55 /im. 

Table [2] Asteroid Pair Photometry 



Asteroid 




11 






9 






;■ 






i 






z 




1986 


JNi 


18.40 


± 


0.03 


16.95 


± 


0.03 


16.46 


± 


0.02 


16.31 


± 


0.01 


16.24 


± 


0.02 


2000 


WXi 67 


21.05 


± 


0.08 


19.65 


± 


0.02 


19.07 


± 


0.02 


18.92 


± 


0.02 


18.82 


± 


0.04 


2001 


MD 30 


18.95 


± 


0.07 


17.62 


± 


0.02 


17.10 


± 


0.01 


16.93 


± 


0.01 


16.89 


± 


0.15 


2000 


NZ 10 


20.15 


± 


0.06 


18.50 


± 


0.02 


17.79 


± 


0.02 


17.61 


± 


0.01 


17.63 


± 


0.03 


2002 


AL 80 


22.51 


± 


0.30 


20.66 


± 


0.03 


20.02 


± 


0.02 


19.83 


± 


0.03 


19.85 


± 


0.07 


1999 


KF 


21.21 


± 


0.12 


19.46 


± 


0.02 


18.77 


± 


0.02 


18.59 


± 


0.02 


18.58 


± 


0.04 


2002 


GP 75 


22.44 


± 


0.27 


21.01 


± 


0.03 


20.36 


± 


0.03 


20.27 


± 


0.03 


20.25 


± 


0.11 


2006 


AL 54 


21.97 


± 


0.20 


20.54 


± 


0.03 


19.88 


± 


0.02 


19.68 


± 


0.02 


19.66 


± 


0.06 


1962 


RD 


17.32 


± 


0.05 


15.54 


± 


0.03 


14.88 


± 


0.03 


14.60 


± 


0.02 


14.63 


± 


0.02 


1997 


CT 16 


21.14 


± 


0.09 


19.37 


± 


0.02 


18.59 


± 


0.02 


18.43 


± 


0.01 


18.44 


± 


0.04 


2000 


RV 55 


22.07 


± 


0.27 


20.31 


± 


0.03 


19.60 


± 


0.02 


19.36 


± 


0.03 


19.41 


± 


0.05 


2004 


RJ294 


23.27 


± 


0.54 


21.77 


± 


0.06 


21.09 


± 


0.05 


20.93 


± 


0.05 


21.00 


± 


0.19 


2003 


sc 7 


22.21 


± 


0.18 


20.60 


± 


0.02 


20.03 


± 


0.02 


19.77 


± 


0.03 


20.15 


± 


0.11 


2000 


GQ113 


20.52 


± 


0.05 


18.89 


± 


0.03 


18.27 


± 


0.02 


18.07 


± 


0.02 


18.24 


± 


0.02 


1983 


WM 


19.59 


± 


0.03 


17.78 


± 


0.01 


17.07 


± 


0.01 


16.87 


± 


0.01 


17.10 


± 


0.02 


2003 


YK39 


22.53 


± 


0.31 


20.56 


± 


0.03 


19.83 


± 


0.03 


19.77 


± 


0.02 


19.81 


± 


0.07 


1999 


TE221 


20.71 


± 


0.08 


19.21 


± 


0.02 


18.66 


± 


0.02 


18.46 


± 


0.02 


18.75 


± 


0.04 


2000 


LU 15 


21.50 


± 


0.11 


19.76 


± 


0.02 


18.99 


± 


0.02 


18.88 


± 


0.02 


19.39 


± 


0.05 


2001 


XH209 


23.30 


± 


0.56 


20.58 


± 


0.04 


19.82 


± 


0.02 


19.57 


± 


0.04 


19.60 


± 


0.10 



Table 2: 

ugriz photometry for 19 asteroid pair members from SDSS DR7 MOC4. The two objects shown in bold constitute 
the only complete pair. 



Table [2] provides multiband photometry shown in Figure [2] f rom the SPSS DR7 MOC4 
for 19 members of 18 pairs from the sub-set of non- family pairs (IPravec et al. 1120091 ). Only 
one pair had observations of both members available in the SDSS data set. 



5. Taxonomic type identification 



While there has been recent and interesting prog ress in automated taxonomic classifica- 
tion from asteroid photometry or spectrometry {e.g. iMarzo et al. Il2009t iMisra &: Busll2008l ; 
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Fig. 2.- 

SDSS five-filter ugriz photometry for 19 objects common to the dynamically young pair 
members and also in the SDSSDR7 MOC 4. The approximate band centers and widths for 
all five filters are shown at the bottom of the Figure. The SDSS filter measurements have 
been connected by lines to simulate full spectra and are normalized to 1.0 at 0.55 fim. 



Bus fc Binzelll2000l ; iHowell et al. 1119941 ) the process is still as much an art as science. In an 
attempt to be consistent (if not rigorous) in our taxonomic classifications we developed and 
employed three closely related methods: 1) taxonomic fitting 2) principal color component 
space location and 3) visual assessment. Our classifications with the three methods were 
correlated with each other and with accepted taxonomic classifications for known objects 
but were not identical. 

The first method selects the best taxonomic 'fit' for an object's photometry to each 
SMASS spectral type j according to the spectral difference 

%=ugriz v °ij ' 
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where Ai and S'jj are the asteroid and SMASS ( iBus and Binzelll2002l ) standard mean flux in 
band % respectively and <ta v <xs y are the errors on the asteroid and SMASS magnitudes in 
the same band. The errors on the mean SMASS spectra range from 0.01 - 0.06 and average 
about 0.03. 

One complication with employing this technique is that the u band center is at 3557 A 
while SMASS spectra extend only over the range 4400-9200 A. Therefore, to compare SMASS 
spectra to the SDSS filter bands in eq. (and in the other two methods described below) we 
performed a linear extrapolation of the SMASS spectra to the u band central wavelength from 
the g and r band centers. We tested several other methods of extrapolating to the u band 
including a quadratic extrapolation to u using g, r and i bands or a quadratic extrapolation 
using all four of g, r, i and z bands. We also tested the result when we simply ignored 
the u band. None of these techniques gave type identifications as good as the simple linear 
extrapolation from g and r even though the lin ear approximation glosses over a possible 



inflection point at 0.42 /xm (jlshiguro et al. 1 120071 ). 



The second method selects the taxonomic type with the smallest distance in principal 
component color space (PC\,PC<2) between the mean for each SMASS type j and the object 
of interest as illustrated in Figure [31 We parametrize the distance Rj simply as 

R) = APCl + APC% (22) 

where APCij and APC 2 j are the difference in PC\ and PC 2 between the object and the 
SMASS class j band average. 

Finally, the third method was a traditional visual assessment of the asteroid's spectrum. 
This time-tested method incorporates elements of both methods discussed above, takes into 
consideration the extrapolation to u band and takes advantage of human perception for 
overall shape matching. In this method, more weight was given to the degree of shape 
matching than to minimizing band center differences. 

In the end, we found that the top two candidate SMASS subclasses from each of the 
three methods were consistent and that visual assessment provided the the most reliable 
taxonomic classification method. 



6. Results and discussion 

6.1. Sub-My cluster member taxonomy 

We observed the brightest member of each of the four sub-My clusters to identify their 
taxonomic type. The spectra or multiband photometry for each object are shown in Fig- 
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Fig. 3.- 

PC2 vs. PC\ for eighteen pair members having SDSS ugriz photometry are shown as circles. 
Filled/open circles lie inside/outside the region generously defining the S-complex. Mean locations 
for SMASS classes are indicated by their taxonomic identification code. We assign Q-class objects 
to the S-complex as discussed in the text. 



ureH (1270) Datura, (16598) 1992 YC2 and (21509) Lucascavin all show classic S-complex 
characteristics in the visible — a 0.75 /im peak and a 1.0 /im absorption band. Datura also 
shows an inflection nea r 0.55 /xm indicating a fresher surface relative to older S asteroids 
with smoother spectra (jlshiguro et al. 1 120071 ; iHiroi et al. 1120061 ) . 



(21509) Lucascavin also shows the 2.0 /im band typical of pyr oxene. Using the techniques 
described in §5] we identify (1270) Datura as an Sk, close to iMothe-Diniz and Nesvorny 



f l2008l ys identification of SI, (21509) Lucascavin also as an Sk, and (16598) 1992 YC2 as 
member of the S-complex. Our visible and IR spectra of (14627) Emilkowalski does not 
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show the 1.0 and 2.0 fxm absorption bands typical of the S-complex and we classify it as 
T type. As our space weathering model only applies to asteroids within the S-complex we 
ignore (14627) Emilkowalski for the remainder of this work. 
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Fig. 4.— 

Visible and near IR spectra of (14627) Emilkowalski and (21509) Lucascavin obtained with 
UH2.2 m/SNIFS and IRTF/SpeX, (1270) Datura visible spectrum obtained with UH2.2 m/SNIFS, 
and SDSS (16598) 1992 YC2 photometry. The spectrum of (14627) Emilkowalski is normalized to 
1.0 at 0.55 //m and the others are offset vertically. The gap at ~ 1.9/xm results from removing a 
sky absorption band. The (1270) Datura visible spectrum, (14627) Emilkowalski and (21509) Lu- 
cascavin IR spectra are all smoothed fits to the data (causing a spurious mismatch to the visible 
spectrum), the others are binned. 



It is unsurprising that (1270) Datura and (21509) Lucascavin were identified in the S- 
complex as both clusters are located in the inner main belt which is dominated by S-complex 
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asteroids. Similarly, (14627) Emilkowalski and (16598) 1992 YC2 are located in the middle 
of the main belt where X and S-complex types are common. 

Having identified three of the sub-My clusters within the S-complex for which we intend 
to examine the effects of space weathering and gardening we obtained data for eight of the 
cluster members as provided in Table |3j 

Table [3j S-complex sub-My cluster members's derived spectral and age data 



Cluster 


Asteroid 


Source 


Slope 


PC X 


Age 








(Reflectance/^m) 




(kyr) 




(1270) Datura 


MDN 


0.355 


0.391 


530 ± 20 




(1270) Datura 


this work 


0.382 


0.414 


!) 


Datura 


(203370) 2001 WY 35 


MDN 


-0.155 


-0.053 


11 




(60151) 1999 UZ 6 


MDN 


0.606 


0.609 


11 




(90265) 2003 CL 5 


MDN 


0.205 


0.260 


11 








mean 


0.305 ± 0.278 






(21509) Lucascavin 


this work 


0.476 


0.496 


300 - 800 


Lucascavin 


(209570) 2004 XL 40 


MDN 


0.314 


0.355 


J) 




(180255) 2003 VM 9 


this work 


0.327 


0.366 


53 








mean 


0.406 ± 0.078 





1992 YC2 (16598) 1992 YC2 MDN 0.083 ± 0.02 0.155 ± 0.027 135 - 220 



Sample Mean 0.36 ± 0.07 511 ± 10 



Table 3: 

Derived color data and age estimates ([Nesvornv et al. 1120061 : iNesvornv and Vokroumickvl 120061 ) for 
members of three S-complex sub-My clusters. PC\ was calculated using eq. [20j Errors on the slope 
and PC\ are not included as the error on the family's color is dominated by the distribution of PC\ 
values within a cluster. The errors on cluster means are standard deviations except for (16598) 
1992 YC2 for which the measurement error is provided since it was the only object observed in the 
cluster. The Sample Mean includes all eight cluster members with the error on the mean. 
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Tableg] Asteroid Pairs 



Asti 


SMASS 


Tholen 


Quality 


Ast 2 


a 


e 


i 


Hi 


H 2 
















(AU) 




(degrees) 






1986 


JNi 


X, Xe 


EMP, EMP 


good 


2001 


XO105 


1.946 


0.0601 


23.710 


13.5 


17.4 


2000 


wx 167 


Xe, T 


EMP, T 


fair 


2007 


uv 


1.909 


0.0613 


23.096 


16.2 


17.1 


2001 


MD 30 


Xe, X 


EMP, EMP 


fair 


2004 


TV14 


1.938 


0.0886 


19.987 


14.9 


17.2 


2000 


NZio 


L, SI 


S, S 


good 


2002 


AL 80 


2.287 


0.1801 


4.097 


14.1 


16.2 


2002 


AL 8 o 


SI, S 


s, s 


good 


2000 


NZ 10 








16.2 


14.1 


1999 


KF 


L, SI 


s, s 


good 


2008 


GR90 


2.327 


0.2339 


1.777 


15.0 


17.2 


2002 


GP 75 


L, S 


s, s 


good 


2001 


UR224 


2.340 


0.1727 


3.865 


15.7 


17.2 


2006 


AL54 


L, SI 


s, s 


good 


2000 


CR49 


2.272 


0.1763 


4.591 


16.8 


14.3 


1962 


RD 


SI, Ld 


s, s 


good 


1999 


RP27 


2.198 


0.1775 


1.129 


13.1 


15.3 


1997 


CT 16 


SI, Sa 


s, s 


good 


2002 


RZ46 


2.186 


0.1672 


4.599 


15.4 


16.4 


2000 


RV55 


SI, Sa 


s, s 


good 


2006 


TE23 


2.657 


0.1026 


2.245 


14.9 


16.8 


2004 


RJ294 


S, Sr 


s, s 


good 


2004 


GH33 


2.268 


0.0981 


4.238 


18.2 


16.7 


2003 


sc 7 


Sk, K 


s, s 


good 


1998 


RB75 


2.264 


0.1114 


7.263 


16.6 


14.6 


2000 


GQ113 


Sq, Sk 


s, s 


good 


2002 


TO 134 


2.324 


0.1319 


5.515 


14.4 


16.3 


1983 


WM 


Sr, Sa 


s, s 


good 


1999 


RC118 


2.320 


0.0790 


5.726 


13.7 


14.6 


2003 


YK39 


Sr, Q 


S, Q 


good 


1998 


FL116 


2.187 


0.0845 


3.736 


18.3 


15.0 


1999 


TE221 


Q, v 


Q, v 


fair 


2001 


HZ32 


2.308 


0.1540 


5.642 


16.5 


15.0 


2000 


LU 15 


V, Q 


V, Q 


good 


1992 


WJ35 


2.313 


0.0701 


5.742 


16.1 


13.7 


2001 


XH209 


A, Sa 


A, S 


good 


2004 


PH 


2.401 


0.2150 


3.638 


15.6 


16.4 



Table 4: Asteroid pair member data for objects appearing in SDSS DR7 MOC4. We de- 
termined SMASS taxonomy for Ast± using ugriz photometry from SDSS DR7 MOC4 as 
described in the text. The first and second ranked SMASS classes (^J5]) are provided along 
with the corresponding Tholen classes and an assessment of that match's quality (degree 
of identification certainty). Semi- major axis (a), ecce ntricity (e), inclination (T) a nd ab - 
solute magnitudes [ H-\,H?) of both asteroids are from IVokrouhlicky and Nesvornyl ( 120081 ) : 



Pravec et al. I ( 120091 ). Members shown in bold constitute the only complete pair. 



6.2. Taxonomy and orbit distribution of asteroid pair members 



None of the members of 36 non-family asteroid pairs ( iPravec et al. I (120091 ) ; IVokrouhlicky and Nesvorny 



( 120081 ). £Q] describes paired aste roids' discovery based on similar orbits.) are available in ei- 
ther the SMASSI or SM ASSII Jl3us and Binzel I2OO2I) spectra databases or the Eight Color 
Asteroid Survey (ECAS) ( lTholenlll989[ ). This is unsurprising considering that the pair mem- 
bers are considerably smaller than the typical asteroid in those s urveys. However, t he 19 
pair members identified in Table H] were found in the SDSS MOC4 (IParker et al. 1120081 ) from 
which we obtained the five-filter solar-corrected ugriz photometry in Table |2j 



Most of the pair members are located in the inner main belt in a region dominated by 
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S-complex asteroids. Their SDSS photometry indicate that they belong to various taxonomic 
types typical of the inner belt including L, S and V classes. Our formal identification of the 
pair member's taxonomy using the methods described in §5] are provided in Table H] which 
shows that we have identified 14 of the 19 pair members with the S-complex in which we 
also include L-class and Q-class. We will examine the colors of these asteroids in the context 
of our space weathering model later in this section. We also identified three X-class, one 
V-class, and one A-class asteroid. 

The taxonomic variety of the pair members is also represented in Figures E] and E] which 
shows that their PC\ color distribution is narrower than the full span of SMASS classes 
while the distribution of their PC2 values extends beyond the SMASS class range. Several 
pair members have PC2 < —0.4. Since PC2 corresponds roughly to a spectrum's curvature 
this indicates an unusually convex spectrum. A couple pair members lie close to the V-class 
region while only one lies in the C-complex. We believe that the dearth of C-complex objects 
is an observational artifact because pair members tend to be small and would be difficult to 
detect with the low albedo of C-complex members in the outer belt. 

The assignation of 1999 TE221 t o the Q class is im portant to our space weathering 



analysis since it has been suggested (IBinzel et al. 1 120041 ) that Q-class objects are actually 



very young, essentially unweathered, S-complex asteroids. Thus, we assume that it is a 
particularly young member of the S-complex with a deep 1/im band as has been predicted 
for young S-complex asteroids. However, in osculating element space it is located close to 
2000 LU15, another member of one of our asteroid pairs from table H] that we assigned to 
the V-class. Both the asteroids lie close to the edge of the Vesta family region as shown in 
fig. [3] (again, in osculating elements). This opens the possibility that 1999 TE221 could be a 
Vestoid with a slightly shallower l/xm absorption band. 

To confirm our Q-class identification for 1999 TE221 and as an additional check on our 
type-identifications we examined whether the pair members are of taxonomic types typical 
of their orbit element phase space region. To do so we identified each pair me mber's five 



nearest orbit element neighbors (using the D-criterion of iNesvorny et al. I (120051) ) in the se t 
of 1175 objects from SMASS that also have osculating orbital elements in AstDyS (j2008 ). 
We found that 17 out of 19 pair members match their nearest neighbor's complex suggesting 
that our identification methods identify the right complex ~ 90% of the time and that the 
pair members are representative of the composition of the main belt region in which they 
are located. 

As mentioned earlier, this test was particularly important for the cases of 1999 TE221 
and 2000 LU15 as both lie on the periphery of the Vesta family region. We identify the 
l st /2 n most likely types for these two objects as Q/V and V/Q respectively. 1999 TE22i's 
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five nearest neighbors include four in the S-complex with one being a Sq and none in the 
V-class. On the other hand, 2000 LU15 has three V-class neighbors. This supports our 
ability to reliably distinguish Q from V. Remember that we place the Q-class within the S- 
complex and, since 1999 TE221 is by far the bluest member of the S-complex pair members, 
its inclusion in our analysis could have a substantial impact on the mean PC\ of the pair 
members and on our measurement of the space weathering rate of S-complex asteroids. The 
effect of including or excluding 1999 TE 22 i in our analysis is described later. 

Our interest in and utilization of the asteroid pairs for the purpose of measuring young 
asteroid surface ages assumes that the pair me mbers are genetically related and fi ssioned 
by some as yet undefined process < 0.5 My ago fjVokrouhlicky and Nesvorny 2008 ). If the 



members of a pair are genetically related asteroids then our expectation is that they will 
display nearly identical spectra. Only one complete pair (2000 NZ 10 and 2002 AL 80 , see 
Table H]) was identified among the 19 pair members available in the SDSS MOC4. Figure [5] 
shows that the colors of the two objects match and therefore supports a genetic origin of the 
pair. Using the taxonomic identification methods of $5] we find that 2000 NZ10 is SMASS 
L-class and 2002 AL 80 is Sl-class - adjacent classes in PC2 vs. PC\ color space as shown in 
Figure |3J Thus, this one line of evidence suggests that asteroid pairs are genetically related. 

Having established the taxonomic composition of the asteroid pairs and their likely 
genetic relationship we would like to examine their taxonomic-orbit distribution — does the 
pair member taxonomy match that of their neighbors in orbit element space? The answer 
to this question could shed light on the relative internal strengths of the different types or 
provide information on the mechanism for asteroid pair creation, i. e. if C-complex asteroids 
split into pairs more frequently it could imply that they are weaker than other types or that 
the pair formation mechanism acts more efficiently on them. Unfortunately, answering this 
question is beyond the scope of this work and we leave it to the future. Instead, we make a 
couple simple observations on the pair's orbit element distribution. 

The axis-inclination structure of the main belt and the pairs is shown in Figure [6] 
revealing that the 36 pairs are distributed in two clumps; a high inclination clump inside 

2.0 AU within the Hungaria family region, a group dynamically protected from perturbations 
by Mars via their high inclination, and a clump on the inner edge of the main belt with 

2.1 AU<a< 2.4 AU. There are also two outliers in the middle belt with semi-major axes 
in the range 2.65 AU< a < 2.75 AU. That most of the pairs are located in the innermost 
main belt is almost certainly an observational selection effect — asteroid pairs are composed 
of small asteroids that are only visible when they are located close to the Earth i.e. in the 
inner main belt. 



The Hungaria clump includes six pairs of which three have SDSS photometry that we 
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Fig. 5.- 

SDSS filter photometry for both members of the only complete dynamical pair in the MOC4. The 
central wavelength for each data point corresponds to the band centers for the ugriz filters except 
for a small horizontal offset for clarity while their width represents the band pass. The data is 
normalized such that a straight line between the g and r data points passes through unity at 
0.55 /im. 



identify as (SMASS) X-complex asteroids — consistent with iGradie and Tesdescol (jl982l )'s 
claim that roughly 70% of the asteroids in the Hungaria region are Tholen E or R class (the 
Tholen E class is contained within the SMASS X-complex). Thus, our identification of three 
SDSS X-complex members in the region is unsurprising and provides further support for our 
taxonomic classification techniques. 
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Fig. 6.— 

Osculating sin(inclination) versus semi-major axis for 36 pair members (large colored points) super- 
imposed on the proper element distribution for main belt asteroids identified in the SDSS MOC4 
(black dots). For semi-ma jor axis < 2.1 AU we show osculating elements for main belt asteroids 
from Astorb (j Astorbl 120081 ) . The 18 distinct pair members with SDSS photometry were identified 
as the following types: violet squares are X-complex, blue triangles are S-complex, green diamonds 
are L-complex, the orange x is V-class, the fuchsia x is Q-class, and the fuchsia + is A-class. The 
18 red asterisks represent pairs for which neither member is present in the SDSS MOC4. The v§ 
and 3:1 resonances are shown for orientation along with the Hungaria family region. 



6.3. Space weathering on sub-My clusters and asteroid pairs 



Figure [7] combines our previous work ( iWillman et al.\ l2008h with the new color-age 
data in this work for S-complex asteroids in three sub-My clusters and eleven sub-My as- 
teroid pairs. The mean PC\ for the three sub-My clusters is 0.36 ± 0.07 — within \a of 
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Willman et al. I (120081 ) 's predicted color of PC\ = 0.31 ± 0.04 for the clusters's mean age of 
440 ± 60 kyr. The good agreement with the prediction lends confidence to the space weath- 
ering model which now agrees with the cluster color-age data over five orders of magnitude 
in age in the decades from 10 5 — 10 10 years. 
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Fig. 7.- 

■PCi color and dyna mically determined ages for S-complex asteroid families adapted from 



Willman et al. I (|2008l ). (The corresponding spectral slope is shown on the right.) The dashed 



curve represents their space weathering model (eq. [3]) extrapolated to the sub-My region. Three 
S-complex sub-My clusters ((1270) Datura, (21509) Lucascavin, (16598) 1992 YC2) are shown in- 
dividually and with their mean value indicated by the red 'Sub-My cluster average' data point. 
The solid curve represents the dual weathering/gardening model fit to the family data including 
the sub-My cluster point. The 'SDSS pairs' point represents the average color of 12 unique S- 
complex sub-My pairs found in SDSS DR7 MOC4. Errors are standard errors on the mean except 
for (16598) 1992 YC2, a single object, for which we provide the measurement error. 
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However, the weighted mean PC\ = 0.49 ± 0.03 color (the error is w eighted error on 
the m ean) for the S-complex young pairs is over 5a redder than predicted by IWillman et al. 
( 120081 ). Excluding the Q-type pair member 1999 TE 22 i discussed in §6.21 increases the mean 
PC\ to 0.54 ±0.04 and the discrepancy to over 7a. T he very young as t eroid pairs clearly do 
not follow the space weathering function proposed by IWillman et al. I ( 120081 ). 



The disparity may be explained in a number of ways: 1) the asteroid pairs do not 
represent a recent breakup of a parent body and are not genetically related or 2) the asteroid 
pairs are the result of a recent breakup but with only partial resurf acing which did not 'res et' 



the space weathering clock or 3) the space weathering model of IWillman et al. I ( 120081 ) is 



either too simplistic or wrong. We examine each of these scenarios in turn: 



1. We consider it unlikely that the asteroid pairs a re not genetically related for two 
reasons. First, the pairs were specifically selected ( IVokrouhlicky and Nesvornyl 120081 ) 
because they are statistically likely to be related asteroids. Furthermore, the colors of 
the pair for which both members exist in the SDSS MOC4 agree extremely well (see 
Fig. [5]). While asteroids that inhabit the same region of the main belt often have 
similar colors the scale of agreement in both the orbit and colors argues persuasively 
for a genetic link between the pair members. 

2. Some pair formation scenarios may not be as violent as the formation of large aster- 
oid families through the catast rophic disruption of a paren t body and may not reset 
the entire surface to zero age. IVokrouhlicky and Nesvornyl (120081 ) cite three possible 
methods of forming pairs : catastrophic colli s ion followed by fragmen t reaccumulation 

YORP induced 
or 



into binary orbits (e.g. iDurda et al. 1 120041 ; iNesvorny et al 



2006) 



Walsh et al. I 120081 ) 



rubble pile spin- up leading to calving of the secondary (e.g. 
YOR P induced angular acceleration of contact binaries leading to their separation 
(e.g. iMerline et al. 1120021 ; iPravec et al. 1 12009c IDurda et al. 1 120041 ) . The first scenario 



will 'reset' the age of the entire surface of both the primary and secondary. The third 
mechanism might leave large portions of the surface undisturbed since the only portion 
that is necessarily exposed is the binary contact region. It is unclear how much of the 
surface would be affected by the second scenario. Therefore there is a clear distinction 
between the possible formation mechanisms of pairs and larger groupings. Asteroid 
families and clusters with more than two members form only by the first catastrophic 
collisional scenario while pairs can also form by the two variations on YORP spin-up. 
This yields at least one pair formation scenario that could leave some of the members' 
surface undisturbed. In this case the globally averaged surface color would correspond 
to a misleading age somewhere between fresh surface and the age of the parent body's 
surface. Taking the color of the pairs at face value and interpreting it in the framework 
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of the IWillman et al. I (120081 ) space weathering model indicates an average surface age 
of ~ 440 My — over 1000 x older than their dynamical age. If we assume that the orig- 
inal parent body's surface was reddened to saturation prior to separation, and taking 
the dynamical ages of < 500 kyr for the pairs at face value such that freshly exposed 
surface is essentially unweathered, then ~ 64% of the surface must be disturbed in 
the pair separation process. Considering that it is unlikely that the surface was fully 
weathered prior to separation allows us to set an upper limit on the fraction of dis- 
turbed surface at < 64%. (Excluding the Q-type object 1999 TE 2 2i only changes the 
upper limit to < 48%.) A gentle bina ry separation d ue to slow YORP spinup may be 
consistent with this scenario. Indeed, iPravec et al. I (120091 ) also provide evidence from 
reconstruction of the initial configuration of the 6070-54827 pair and rotation rate ob- 
servations that the non-family pair formation process is a gentle event. If we envision a 
bi-lobed asteroid gradually accelerating in angular velocity and finally fissioning at the 
neck that joins the lobes then it seems reasonable that the portion of dist urbed surface 



would be < 64%. In the binary/pair formation mechanism proposed by IWalsh et al. 



( 120081 . 120061 ) an asteroid's polar surface material migrates to the equator as the ob- 
ject's rotation rate increases and eventually flies into orbit around the parent body 
where it reaccumulates into a satellite. The primary and satellite eventually separate 
and their orbits evolve dynamically. Our impression is that this model would generate 
fresh (blue) surface on both the primary and satellite in conflict with our observation 
of reddish surfaces on the asteroid pairs. However, it is not difficult to envision a slow 
migration process that allows material to weather on the primary's surface before being 
shed and accumulating into a secondary object. 



3. It is possible that t he space weathering mo del o f IWillman et al. I (120081) that built upon 
the earlier work of iNesvorny et al. I ( 120051 ) and iJedicke et al. I (J2004J) is simply wrong; 
that the apparent change in color of S-complex asteroids with age is a statistical fluke 
or due to some other underlying effect (though obvio u s pos sibilities were considered 
in detail in the early works). However, iParker et al. I ( 120081 ) confirm the weathering 
effect in an independent updated analysis of the SDSS DR7 MOC4 data. We consider 
it more likely that the space weathering mechanism is more complicated than simply 
affecting the average spectral slope (or PC\) as a function of time. For instance, it is 
well known that space weathering affects not only the slope of the spectrum but also 
the depths of the 1 /im and ultraviolet absorption bands and the surface albedo. We 
consider it not only possible but likely that these effects occur during space weathering 
at different rates. The apparent redness of the young asteroid pair members relative 
to the expectation of the simplistic space weathering model could indicate that a 'fast' 
space weathering process takes place in < 10 5 years. However, this scenario requires 
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a rather contrived sequence of events. Consider the three color change processes: 
decreasing depth of the ultraviolet band shortward of 0.4 /im, decreasing depth of the 
1 /im absorption band and continuum reddening between these two bands. The first 
process is the only one that produces a bluer color. Therefore, accounting for the 
anomalous redness of the pair members requires the unlikely scenario that one of the 
latter two processes dominate on short time scales which is then belatedly overwhelmed 
by the first process which is then finally dominated by the third. 



Given the disagreement between the asteroid pair color-age and the IWillman et al. I (120081 ) 
space weathering model, and considering our enumerated arguments above, we continue with 
our analysis under the assumption that we can ignore the colors of the asteroid pairs in this 
new determination of the space weathering and gardening rates. 

First, we fit all the S-complex color data including the sub-Myr clusters but not the 
asteroid pairs to the 'old' space weathering fun ction of eq. [3J Considering the good agreement 



between the predicted (IWillman et al. 1 120081 ) and observed colors of the sub-Myr cluster 



members it is unsurprising that the fit including the new data matches the previous fit 
in all four parameters to within 1-a with PCi(O) = 0.34 ± 0.02, APCj = 0.28 ± 0.05 , 



t = 700 ± 270 My, a = 0.58 ± 0.17. However, as we observed in IWillman et al. I (120081 ). 
fitting the color-age data to the form of eq. [3] suffers from multiple and wide minima in the 
fit-parameter space. Furthermore, the function does not explicitly separate the weathering 
and gardening effects. 

On the other hand, we found that fitting the same color-age data to our new function 
that incorporates both space weathering and gardening (eq. [HD is better behaved because the 
solution space does not show multiple local x 2 minima. The best fit (lowest x 2 ) including 
the sub-My clusters (but, again, not the asteroid pair colors) yields PCi(0) = 0.37 ± 0.01, 
APCi = 0.33 ± 0.06, t w = 960 ± 160 My, r g = 2000 ± 290 My as shown by the solid curve 
in Figure [3 

Note that it is not correct to compare the new t w or r g to the old r value because of the 
complicating and non-physical use of a in the generalizing exponent in the old functional 
form. In essense, the 'old' value was an effective weathering time that combined the effects 
of both regolith gardening and weathering. Since the effective weathering time can be shown 
to be equivalent to the time corresponding to the inflection point on the new weathering 
function, Figure [7] shows that the two models are in good agreement. 

The old and new APCi are not strictly comparable either. Formerly in the weathering 
only case the entire surface would eventually reach the color PC\ = PCq + APC. The 
weathering/ gardening case will produce a lower equilibrium value with PC\ < PCq + APC. 
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The new space weathering-only time frame of ~1 Gy is consistent with the 'slow weather- 
ing' measured in lab-based measurements as discussed in the introductions. Th e 'slow' school 
includes our result of ~1 Gy based on space observations of S-complex families, ISasaki et al 



equivalent value of 700 My at 2.6 AU based on pulsed-laser irradiated silicate pellets, 
and lPieters et al. I ( 120001 ) 's estimate of 100 — 800 My for lunar surfaces based on crater counts 
and spectra and color-terrain correlations on (243) Ida (jVeverka et al. II1996I ). But these slow 
space weathering results di sagree with the ' fast w eathering' sch ool that includes lab- based 



ion irradiation results from IVernazza et al. 
5 kyr 



((20091) ( < 1 My) and lLoeffler et al. I fl2009D's (~ 
and iTakatol (120081 ) 's <450 kyr time scale based on a shallow 1 /im absorption band 
observed on (1270) Datura. We do not have an explanation for the discrepancy between the 
fast and slow weathering results. 



6.4. Gardening Time 



In §3] we calculated the gardening time for asteroid regolith as a function of target 
diameter as shown in Figure. [TJ Over a broad range of asteroid sizes the gardening 
time scale is 200 — 300 My years. T his result is in dramatic disagreement with a simi- 



lar calculation by iMelita et al. I (120091 ) that was based on a collisional cratering model by 
Gil-Hutton and Lazzarol (120021 ) and yielded a time scale of 2000 years for resurfacing Trojan 
asteroids. It is difficult to reconcile the two results that differ by five orders of magnitude. 
Two factors that we are aware of that explain some of the discrepancy are 1) that they used 
a Trojan impact probability double that of the main belt and 2) the slope of the Trojan 
size-frequency distribution was assum ed to be —3 whereas the data we used had slope of 



~ —2.20 appropriate to the main belt (IBottke et al. Il2005bl ). But these two differences only 
account for a fraction of the difference between our results. 

Figure [1] also shows that our measured gardening rate from the S-complex asteroid 
color-age relationship is an order of magnitude different from our calculated rate based on 
impacts. Is it possible to reconcile this difference? 

There is considerable uncertainty in the various terms involved in calculating the impact 
gardening time from eq. [19] but we have done our best to select the best contemporary 
values in each case. While the gray region on the figure represents the formal 1-er error on 
the calculation based on the reported errors on each input quantity there is considerable 
unreported systematic error in the calculation associated with its sensitivity to the input 
parameters and functions. In particular, we examined the sensitivity of the calculation to 
the: 
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asteroid size frequency distribution, N(D) 

specific shattering ene rgy fun c tion, Q* s , that determines the largest non- disruptive 



impactor (we used the iMeloshl ( 119971 ) function as a comparison), 

• impact speed and, 

• D min , the smallest impactor that creates ejecta 

In each case we varied the input parameter or function over a range of of 2 to 4 in 
each direction. In most cases the gardening rate for large target objects is only slightly 
affected. The rate for objects ~1 km diameter changed by a factor of about two or slightly 
more. The two most important factors in determining the gardening rate are 1) D m i n and 2) 
the amount of area covered by crater ejecta. Increasing D min or decreasing the ejecta area 
coverage both work to increase the gardening timescale. Leaving all other parameters at 
their nominal values, increasing D min to ~30 m is sufficient to increase the gardening time 
to our measured value of ~2 Gy. A similar result is obtained by decreasing the diameter of 
the ejecta field by a factor of two (the affected area by a factor of four). 



In a related observation, I Chapman! (120051 ) and others note the unexpected paucity of 



craters <200 m diameter on (433) Eros. T aken at face va lue the observation would imply a 



deficit of impactors < 20 m in diameter but lO'Brienl ( 120091 ) has shown that no such deficit can 



exist because it would in turn generate a 'wave' in the observed size-frequency distribution 
for larger asteroids that is not observed. 

Thus, we have identified two independent problems — the mismatched measured color- 
age and impact-calculated gardening times and the crater deficiency on (43 3) Eros — that can 



both be reso lved if impactors < 20 m do not leave a crater record. While iRichardson et al. 



( 120041 . 120051 ) have proposed that seismic shaking can erase the small craters the mechanism 
can only explain the crater shortage not the gardening time mismatch. This is because 
any process that increases the gardening rate (i.e. seismic shaking) shortens the calculated 
impact gardening time and worsens the mismatch with measured color-age gardening rate. 

If we believe that there is no deficit of impactors of < 20 m diameter and that there is a 
deficit of craters and regolith gardening caused by impactors in the same size range then it 
must be the case that the impacts take place without creating craters in the 'normal' manner. 
Perhaps the surfaces of small asteroids 'absorb' the small impactors in an inelastic collision 
without creating a crater an order of magnitude larger or ejecting significant material. Like 
throwing a small stone into a bee's empty honeycomb. 
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7. Summary 

We have combined data from five aster oids obtained on Keckll/ESI, UH 2.2/SNIFS 



and fRTF/SpeX, spectra of six asteroids from iMothe-Diniz and Nesvornyl ( 120081 ) and SDSS 



archive photometry on 19 asteroid pair members to investigate the earliest stages of space 
weathering on asteroids aged less than one million years, including sub-My old asteroid 
clusters and asteroid pairs. 



Willman et al. I (120081 ) predicted that the color of these young S-complex asteroid's 
surfaces would be essentially identical to freshly exposed regolith with PC\ = 0.31 ± 0.04. 
Our measured color data for asteroid clusters with ages in the range 10 5 — 10 6 years agrees 
with their prediction with a mean PC\ = 0.36 ± 0.07. 

For the first time, we extended the space weathering model to explicitly include the 
effects of regolith gardening which restores weathered S-complex asteroid surface to its orig- 
inal blue color. Fitting the refined model to the former data along with the new sub-My 
asteroid cluster data point gives separate characteristic times for weathering and gardening 
of t w = 960 ±160 My and r g = 2000 ± 290 My respectively. The new results suggest that 
fresh S-complex asteroids have a blue surface color with PCi(0) = 0.37 ± 0.01 and would 
redden by APC\ = 0.33 ± 0.06 over a very long period of time in the absence of gardening. 
The presence of gardening produces an ultimate equilibrium color of PC^(o o) = 0.59 ±0.06. 



The n ew weathering time scale and colors are consistent with the results of IWillman et al 



( I2008f ) 



Our new data for asteroid pair members with ages < 10 5 years is redder than predicted 
by more than 5a. Assuming that the space weathering model is correct, the discrepancy 
could be due to the pair production formation mechanism — if the asteroids's surfaces are 
not fully recoated during a gentle separation then the surface age for these pairs could be 
much older than the dynamical age since separation. This explanation requires that < 64% 
of the asteroids's surface is disturbed during the pair formation event providing an interesting 
test in comparison to simulations of pair formation. An alternative explanation that there 
are 'fast' and 'slow' weathering processes occuring on S-complex asteroid surfaces does not 
seem plausible. 

We independently calculated the gardening rate on main belt asteroids from basic prin- 
ciples including the affects of the asteroid size-frequency distribution, impact rates, crater 
and ejecta formation. The calculated gardening time scale of ~ 270 My for a D ~3 km 
asteroid typical of those in our sample is in stark contrast to the color-age measured value of 
~ 2000 My. Two scenarios that can reconcile the calculated and measured gardening rates 
are 1) impacting asteroids of < 20 m diameter do not produce craters in the standard manner 
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and 2) the ejecta field is much smaller than standard models would suggest. Perhaps the 
smaller asteroids are absorbed in inelastic collisions with a 'honeycomb'-like surface. This 
mechanism also provides an explanation for the paucity of small craters on (433) Eros. 
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